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We explore the nature of glass-formation in variable spatial dimensionality (d) based on the 
generalized entropy theory, a synthesis of the Adam-Gibbs model with direct computation of the 
configurational entropy of polymer fluids using an established statistical mechanical model. We find 
that structural relaxation in the fluid state asymptotically becomes Arrhenius in the d ^ oo limit 
and that the fluid transforms upon sufficient cooling above a critical dimension near d = 8 into a 
dense amorphous state with a finite positive residual configurational entropy. Direct computations 
of the isothermal compressibility and thermal expansion coefficient, taken to be physical measures 
of packing frustration, demonstrate that these fluid properties strongly correlate with the fragility 
of glass-formation. 

PACS numbers: 64.70.Q- 


Growing evidence indicates that molecular packing 
plays an essential role in both the physics of crystalliza¬ 
tion [IHi and glass-formation prompting a con¬ 

sideration of variable spatial dimensionality (d) as a con¬ 
ceptual probe of these solidihcation processes. In par¬ 
ticular, the “decorrelation principle”, recently proposed 
by Torquato and Stillinger [ll, [T^, states that uncon¬ 
strained correlations in hard hyperspheres diminish with 
increasing d and vanish in the d —> oo limit. Likewise, 
recent simulations have indicated that certain of¬ 

ten cited aspects of glass-formation, such as the dynamic 
heterogeneity and the decoupling between structural re¬ 
laxation and diffusion, become diminished at elevated d. 
While prior studies |lll - [l^ have focused mainly on how 
critical packing fractions associated with glass-formation 
in hard hyperspheres vary with d, little is known concern¬ 
ing the d dependence of the characteristic temperatures 
and fragility of glass-formation in molecular fluids. No 
prior work has considered complex fluids, such as poly¬ 
meric liquids, where many molecular parameters may be 
tuned to control the nature of glass-formation. 

In this Letter, we consider the glass-formation of a 
model polymer glass-forming (GF) liquid in d dimen¬ 
sions. Our approach is based on the generalized entropy 
theory (GET) [2^, which is a combination of the lat¬ 
tice cluster theory (LGT) [2l| and the Adam-Gibbs (AG) 
theory linking the configurational entropy Sc (i-C., 
the entropy devoid of the vibrational component) to the 
structural relaxation time Tq,. The LOT employs a d- 
dimensional hypercubic lattice model, whose use facili¬ 
tates the development of an expansion about the mean- 
field limit of infinite d. The LGT thus naturally provides 
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a framework for considering the d dependence of poly¬ 
mer thermodynamic properties, thereby making the GET 
suitable for exploring glass-formation in d dimensions. 

As a primary result, we find that structural relaxation 
in the fluid state asymptotically becomes Arrhenius in 
the d oo limit [23|, a limiting behavior that is a 
dynamical counterpart of the decorrelation principle of 
Torquato and Stillinger Bin. Moreover, our theory 
indicates that the fluid transforms upon cooling above 
a critical dimension near d = 8 into a dense amorphous 
state with a finite positive residual configurational en¬ 
tropy. We show that the “kinetic fragility”, determined 
from the temperature (T) dependence of Tq,, decreases 
with d, while the “thermodynamic fragility”, defined as 
the relative rate of change of Sc with respect to T, in¬ 
creases with d, in accord with recent simulations for soft- 
sphere liquids for d = 2, 3, and 4 0. Analyses of the 
isothermal compressibility and thermal expansion coeffi¬ 
cient, taken to be physical metrics of packing frustration, 
indicate that molecular packing becomes more efficient 
in higher d, and, as anticipated, these thermodynamic 
properties strongly correlate with the d dependence of 
the fragility. We thus confirm the qualitative geometrical 
origin of the reduced kinetic fragility upon increasing d. 


We consider a d-dimensional melt of chains with the 
structure of polypropylene (PP) B. For comparative 
purposes, a common parameter set is used in all d: the 
microscopic cohesive energy is e = 200 K, the bending 
energy is Ef, = 800 K, and the polymerization index is 
Nc = 8000. By employing these parameters, the GET 
calculations for constant pressure (P = 0 MPa) polymer 
melts for d = 3 yield the glass transition temperature 
Tg and fragility which are typical of synthetic GE 
polymers. 

The temperature dependence of the configurational en¬ 
tropy is central to the AG model. Hence, we first focus 
on this quantity in d dimensions. The GET considers 
the configurational entropy density Sc (i-C., the configu- 
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FIG. 1. Main: Sc/ks as a function of T for various d at 
P = 0 MPa. Dotted lines highlight where Sc is negative, an 
unphysical behavior that is probably an artifact of high T 
expansion in the LCT for low d and a mathematical feature 
that disappears for higher d. Dashed lines denote the low T 
constants of Scjks for each d. Inset: s^/sc as a function of 
(Ta — T)ITa for various d. 


rational entropy per lattice site HI) at constant P, 
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where F is the Helmholtz free energy 
is the total volume of the system with Ni and I^eii the 
number of total lattice sites and the volume associated 
with a single lattice site, and Np is the number of poly¬ 
mer chains. While Well in d dimensions may be defined in 
terms of the volume Well = Oceii ^ d-dimensional cube 
with Cceii the edge length, basic quantities, such as Boltz¬ 
mann’s constant ks, cannot be determined in d > 3. We 
surmount this problem by performing all calculations at 
P = 0 MPa in all d since this condition completely elim¬ 
inates problems in estimating Well and in d > 3. 

Figure 1 displays the calculated ScjkB as a function of 
T for various d. As in d = 3, Sc for larger d displays a 
shallow maximum s* (crosses in Fig. 1) at a high T that 
defines the onset temperature Ta of glass-formation in 
the GET that signals the onset of non-Arrhenius behav¬ 
ior of Ta [1^, and the product ScT exhibits an inflection 
point (asterisks in Fig. 1) at a lower T that uniquely 
defines the crossover temperature W in the GET that 
separates two regimes of T with qualitatively different 
dependences of Ta on T [2^. Evidently, Sc formally be¬ 
comes negative (dotted lines in Fig. 1) at low T and 
low d in our mean-field theory, an unphysical behavior 
that is probably an artifact of high T expansion im¬ 
plicit to the LGT [2^. However, taken at face value, 
the point at which Sc vanishes signals a divergence of Ta 
in the AG model. The existence of an equilibrium solid¬ 
ification transition in a disordered material implies that 
the transition must have the mathematical character of 


a second-order phase transition, exactly the interpreta¬ 
tion made long ago for the glass transition by Gibbs and 
DiMarzio [^ . 

The GET generally predicts the presence of a low T 
plateau in Sc (denoted by Scy, see the dashed lines in 
Fig. 1). This plateau increases linearly with ln(d/2) 
and becomes positive beyond a critical dimension dc [27j| , 
and hence, an “ideal glass transition temperature” Tq, 
where Sc formally vanishes, only exists for d < dc in the 
GET. We further find that dc is independent of e and 
Eb but varies with Nc for a fixed molecular structure 
and that dc is near d = 8 for the PP structure ■ Some 
recent models of glass-formation [il[i3 predict a critical 
dimension of dc = 8 above which the Stokes-Einstein 
relation holds, a prediction that might be related to the 
behavior displayed in Fig. 1. 

The evidence for a nonzero residual entropy in the 
glassy state has been a matter of long-standing discus¬ 
sion . According to the AG theory, the appearance of 
a positive constant Sc implies that does not diverge at 
T > 0 K and that structural relaxation becomes Arrhe¬ 
nius again at low T, albeit with a much higher activation 
free energy than that above Ta- The ratio of the Ar¬ 
rhenius activation free energies in the low and high T 
Arrhenius regimes above dc can be calculated as s*/sc,r, 
a quantity that is evidently larger than unity and that 
decreases with d as s*/sca = 1/(A —H/ln(d/2)] with the 
fitted constants, A = 0.325 and B = 0.505 Since 
Ta becomes astronomically large in the low T Arrhenius 
regime, the thermodynamic state in this regime can be 
considered to be a “glass”. Notably, Sc appears to ex¬ 
hibit a step drop towards Sc,r near W for large d, sug¬ 
gesting that the fluid transforms into a glass at a finite T 
upon cooling through a first-order phase transition in the 
d —oo limit. Interestingly, the “tiling model” of g lass- 
formation, developed by Weber and Stillinger 1^ . 
predicts the glass transition to be first-order. 

We now consider the T dependence of Tq, in d dimen¬ 
sions, which is obtained in the GET through the AG 
relation [^ . 

Ta = Too exp[eAG/ (ksT)], Sag = A^(sc/sc), (2) 

where Too is the high T limit of Ta, eag is the AG ac¬ 
tivation free energy, and A/r is the T independent acti¬ 
vation free energy at high T above Ta where the Arrhe¬ 
nius dependence of Ta holds. Too in d = 3 is taken in 
the GET as a “default” value of 10“^^ s. Since Too is 
expected to be non-universal in variable d, we consider 
only the dimensionless relaxation time = Ta /tqo . The 
GET estimates Au, in d = 3 from the empirical relation, 
Ap/ks = 6 Tc [i^- Simulations for d = 2, 3, and 4 BUI 
suggest that Ap increases approximately in proportion 
to d. This phenomenology is quite understandable from 
classical arguments by Eyring relating Ap to the co¬ 
hesive interaction strength of the fluid, which scales in 
a hypercubic lattice model in proportion to the coordi¬ 
nation number (2d). The scaling of Ap with 2d implies 
that Ap in the GET obeys Ap/ks = 2d Tc^d= 3 , where 
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FIG. 2. log(rr) as a function of T for various d. 



FIG. 4. Main: ScTjkB as a function of TjTx for various d. 
Inset: Kt as a function of d. 


Tc,d =3 is the Tc in d = 3. 

Before analyzing the T dependence of we investi¬ 
gate the T dependence of s*/sc in d dimensions, which 
describes the multiplicative increase of £ag- The inset 
to Fig. 1 displays s*/sc for various d as a function of 
{Ta — T)/Ta-, a reduced temperature that measures the 
degree of quench into the “conjested” regime of glassy 
dynamics. The growth of sag with increasing d and, cor¬ 
respondingly with the scale of collective motion, weakens 
for fixed {Ta — T)/Ta- Evidently, the GET predicts that 
structural relaxation in the fluid state becomes Arrhenius 
in the d —> oo limit where s * /S r approaches unity asymp¬ 
totically. Previous work indicates that s*^/sc for 
d = 3 exhibits a parabolic dependence on {Ta —T)/Ta, 

s*Js, = 1 + Cs[{Ta - T)ITaY, T,<T< Ta. (3) 

Cs measures the extent of cooperative motion at high T 
in the GET [2^, and we find that Eq. (3) holds well in 
{Ta — 50, Ta) in d > 3. 



FIG. 3. Avft, Tvft/Ta, and Cs as a function of d. 


Figure 2 displays the T dependence of log(rr.) for vari¬ 
ous d, along with the fits to the Vogel-Fulcher-Tammann 
(VFT) equation, 


Tr = Tq exp 


1 

KypT{T/TypT — 1 ). 


(4) 


where tq, ATvft, and Tvft designate the high T limit 
of Tj., the “kinetic fragility” parameter quantifying the 
degree to which relaxation is non-Arrhenius, and the 
VFT temperature where Tr extrapolates to infinity, re¬ 
spectively. Previous analyses for d = 3 [2^ have shown 
that the VFT equation describes Tq calculated from the 
GET quite well within the T range between T^ and Tg. 
Since we are unable to define Tg in d > 3 in terms of a 
characteristic timescale, VFT fits are performed for the 
T range between 2k and a lower temperature Ti at which 
log[Tr(T;)] = 63 Q. The d dependence of ATvft is de¬ 
picted in Fig. 3, along with Cs and the ratio Typj /T a, 
a measure of the “breadth” of glass-formation [^. All 
these three metrics for glassy dynamics decrease with d. 

Following the simulation studies 0, we additionally 
explore how d alters the “thermodynamic fragility” pa¬ 
rameter Kt, determined from the slope of the T depen¬ 
dence of the product ScT, 


ScT/kB = Kt{T/Tk - 1), (5) 

where Tr- is the Kauzmann temperature, estimated by 
linearly extrapolating ScT to zero over the same T range 
where KypT is obtained (see Fig. 4). We prefer to 
call Kt the “low T thermodynamic steepness” param¬ 
eter because this parameter bears no direct relation to 
the strength of the T dependence of Tq,. The inset to 
Fig. 4 indicates that Kt increases with d, in agreement 
with the trend found in recent simulations of soft-sphere 
GF liquids in d = 2, 3, and 4 [^. Reference Q also shows, 
based on limited data from simulations, that KypT seems 
to weakly increase with d in the same model |34l | . 
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FIG. 5. Correlations between TFvft, Tvft/Ta, Cs and (a) 
Kr{TA) and (b) ar{TA). 


Our theory, thus, predicts that i^vFT varies in the op¬ 
posite direction as Kt when d grows. An opposite de¬ 
pendence of ATvft and Kt is also found in 3-dimensional 
simulations of GF liquids in which particles interact with 
a modified Lennard-Jones potential (35j |. where the po¬ 
tential softness is tuned to alter the fragility and where 
the inverted variations of Kt and iFvFT have been at¬ 
tributed to variations in A/i. Apparently, the same mech¬ 
anism also explains our results since Ap, scales with 2d 
so that Ap, can be the dominant contribution to i^vFT 
for high d. 

We now seek what physical factors govern the fragility 
of glass-formation. Building on the earlier GET calcu¬ 
lations for d = 3 indicating that the fragility of poly¬ 
mer GF liquids primarily arises from packing frustra¬ 
tion [ 2 ^, we explore physical metrics that might best 
describe packing frustration. We consider the isother¬ 
mal compressibility kt = —(1/F) {dV/dP)\j, and the 
thermal expansion coefficient ap = (1/F) {dV/dT)\p as 
possible physical measures of packing frustration. In par¬ 
ticular, we analyze Kr = KpikHT)/ p (where p = p Well 
with tp the polymer filling fraction [2l|) and = Tap, 
since these reduced quantities are found to strongly cor¬ 
relate with fragility in d = 3. (A more detailed discussion 
for d = 3 will be presented elsewhere.) As expected, as 
d grows, the T dependence of Hr and weakens and 
their overall magnitudes greatly decreases, thereby re¬ 


flecting the diminished packing frustration with increas¬ 
ing d [^. As an illustration. Fig. 5 exhibits the presence 
of strong correlations between various dynamic metrics 
(ATvft, Typt/Ta, and Cs) and Hr{TA) and ar{TA), im¬ 
plying that the d dependence of fragility can be predicted 
from a knowledge of the physical measures of packing 
frustration. Therefore, both the reduced fragility with 
increasing d and the corresponding return to Arrhenius 
relaxation are consequences of the reduction of packing 
frustration in higher d where space is more “open” and 
packing constraints are more weakly “felt”. 

In summary, our exploration of polymer glass- 
formation in variable dimension indicates some clear 
trends as the dimensionality is varied. Perhaps most basi¬ 
cally, we find that structural relaxation in the fluid state 
asymptotically becomes Arrhenius in the d —>■ 00 limit 
and that the fluid transforms upon cooling above a crit¬ 
ical dimension near d = 8 into an amorphous state with 
a finite positive residual configurational entropy. The 
transition from a fragile glass-forming liquid in low di¬ 
mensions to an ideally strong liquid in the d —>■ 00 limit 
indicates that the kinetic fragility decreases with increas¬ 
ing dimension. In accord with simulations Q , our theory 
predicts that the thermodynamic fragility Kp grows at 
elevated d. Kp fails to be a fragility parameter, and 
we suggest that Kp be termed the “low temperature 
thermodynamic steepness” parameter. By establishing 
strong correlations between the isothermal compressibil¬ 
ity and thermal expansion coefficient and the fragility, we 
suggest that these thermodynamic properties should pro¬ 
vide physical metrics of packing frustration for predicting 
fragility variations in glass-forming liquids. 
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Dimension 


SI. MEAN-FIELD ESTIMATION OF THE 
CONFIGURATIONAL ENTROPY AT LOW 
TEMPERATURES AND DEPENDENCE OF THE 
CRITICAL DIMENSIONALITY ON MOLAR 
MASS 

The residual configurational entropy density Sc,r [i-e., 
the low temperature (T) constant of the configurational 
entropy density Sc] can be obtained by analyzing the ther¬ 
modynamics of the system in the limit where all lattice 
sites are occupied by chain segments. Hence, we focus 
on the low T regime of glass-formation where the poly¬ 
mer filling fraction ip (defined as the ratio of the total 
number of chain segments to the total number of lattice 
sites dH, approaches unity in order to estimate how 
the critical dimensionality, above which Sc ceases to van¬ 
ish, depends on molecular parameters. For low T (e.g., 
T = 10 K), the factor exp(—/JF!;,) [where /? = l/(fc_B'r) 
and Eb is the bending energy] in the LCT free energy 
nearly vanishes for large Eb [e.g., the present work uses 
Eb = 800 K, leading to exp(—= exp(—80) at T = 10 
K]. In this limit, all terms multiplied by exp(—in the 
LCT free energy disappear, and then taking ip = 1 leads 
to the greatly simplified expression for the free energy. 



where M is the molar mass (defined as the total number 
of chain segments in a single chain), L designates the 
number of sub-chains in a single chain [2lL l36l | , and Ct 
has the following form [s^: 

G = G,o + G,,(/3e) + G.,2(/3e)^ (S2) 

where e is the nearest-neighbor van der Waals interaction 
energy. However, Ci,e^ and are now indepen¬ 
dent of r, and, hence, 

= + (S3) 

Consequently, the residual configurational entropy den¬ 
sity, Sc,r = - dfr/dT\^, emerges as 

^=-—1 ^ 
ks M ^ Mj 

i—1 ^ 

Substituting Eq. (S3) into Eq. (S4) leads to the relation, 
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+ ^[G,o-G,,2(/3e)2] . (S5) 

i=l 



FIG. SI. Residual configurational entropy density Sc.r/fcs as 
a function of d for various polymerization indices Nc for a melt 
of PP chains. The solid lines are fits to the equation SrjkB = 
Aln(d/2) -I- B with A and B being the fitting parameters. 
The critical dimensionality dc, where Sc,r starts to be positive, 
increases with Nc and saturates at d = 10 for sufficiently large 
Nc 

Taking advantage of the fact that X]i=i = 0 [^ . 
Sc,r finally appears in the simple form: 



Since Ci^ is now a function of c? (2 = 2d) and a 
set of geometrical indices that are determined from the 
monomer structure and molar mass M, Sc,r is indepen¬ 
dent of T, e, and Eb- Figure SI displays the d dependence 
of Sc,r for various polymerization indices Nc for a melt 
of PP chains. The “critical dimensionality” d^ where 
Sc,r starts to be positive, is found to increase with Nc, 
saturating at c? = 10 for large Nc- The d dependence of 
Sc,r/kB in each case can be well described by the equa¬ 
tion, 

Sc,r/kB = Ahi{d/2) + B, (S7) 

where A and B are the fitting parameters. Figure SI il¬ 
lustrates this variation. The filling fraction at the point 
at which Sc formally vanishes is somewhat less than unity, 
but this simple analytic argument is sufficient for esti¬ 
mating the Nc dependence of the critical dimensionality 
above which Sc no longer vanishes at low T. 

S2. DIMENSIONAL DEPENDENCE OF THE 
RATIO OF THE ACTIVATION FREE ENERGIES 
IN THE HIGH AND LOW TEMPERATURE 
ARRHENIUS REGIMES 

As discussed in the main text, the ratio s^/sc^r de¬ 
fines the relative magnitude of the activation free energy 








S2 



d 

FIG. S2. Activation free energy ratio sl/sc,r as a function 
of d for d > dc- The line is a fit to the equation sl/sc,r = 
1/[A — i3/ln(d/2)] with the fitted constants, A = 0.325 and 
B = 0.505. 

for relaxation in the high and low T Arrhenius regimes 
above dc that we separately characterize as being simple 
fluid-like and glass-like because of the large difference of 
relaxation times in these different regimes of T. The GET 
allows us to calculate this activation free energy ratio as 
a function of d, and we illustrate our findings in Fig. S2. 
Curiously, the activation free energy ratio falls with in¬ 
creasing d, indicating that the dynamics of the high and 
low T regimes actually becomes similar in high d, where 
the activation free energy ratio saturates at a minimal 
value of about 3. The d variation of the activation free 
energy ratio is well represented by the approximant: 

s*/ Sc,r = ^/[A- B/\Ti{d/2)], d>dc, (S8) 

where the fitting parameters are estimated to be A = 
0.325 and B = 0.505. Because the activation free energy 
in the high T regime is large for high d, increasing this 
quantity by a factor of 3 in high d translates into an 
extremely large increase in relaxation times in the low T 
Arrhenius regime. 

S3. DIMENSIONAL DEPENDENCE OE THE 
REDUCED ISOTHERMAL COMPRESSIBILITY 

AND THERMAL EXPANSION COEFFICIENT 

Materials that pack more efficiently can be expected 
to have a lower and more weakly T dependent isother¬ 
mal compressibility (kt) or thermal expansion coefficient 
(ap), while those exhibiting a high degree of packing frus¬ 
tration due to the complexity of the molecular structure 
should exhibit a relatively large and relatively strongly 
T dependent kt or ap. Previous analyses of GET 
calculations in d = 3 reveal the existence of a very strong 
correlation between the filling fraction at Tg and the 
fragility when the molar mass, cohesive interaction, or 



T/Ta 



FIG. S3, (a) Main: Reduced isothermal compressibility Kt = 
KriksT) / p as a function of TITa for various d. Grosses, 
asterisks, and circles indicate the points corresponding to Ta, 
Tc, and Tvft, respectively. Data are presented in the T range 
(Tvft, Ta) of interest in the present work. Inset: Same plot 
but shown in a linear scale representation, (b) Main: Kr{TA) 
as a function of d. Inset: Kr{Tc) and Kr(2vFT) as a function 
of d. Lines in the set indicate the exponential fits, nr{Tx) = 
Ax exp{—Bxd) {x = c or VFT) with Ax and Bx the fitting 
parameters. We obtained the fitting estimates: Ac = 11.7 
and Be = 1.05 and Avft = 16.3 and Bvft = 1.77. 


pressure is varied. Moreover, molecules in fragile poly¬ 
mer fluids are predicted to pack less efficiently in the 
glassy state than those in stronger glass-forming poly¬ 
mer fluids. These observations deduced from the GET 
calculations indicate that packing efficiency is a central 
factor for understanding the fragility of glass-formation. 
This line of reasoning naturally leads to the expectation 
that Kt and ap correlate strongly with the fragility. 

Because our extensive analyses in d = 3 reveal 
that the dimensionless isothermal compressibility, = 
KTikpT)/p with p = (p 14eU, and the dimensionless ther¬ 
mal expansion coefficient, ar = Tap, strongly correlate 
with the fragility of glass-formation, we consider these 
reduced quantities in d dimensions. Figures S3 and S4 
indicate the GET model predictions for k^- and ar as a 
function of T/Ta for various d. Evidently, the strength 
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of the T dependence of both and weakens with in¬ 
creasing d, and both Kr and ar progressively decreases 
with d for fixed reduced temperatures T/Ta- The insets 
to Figs. S3 and S4 indicate that Kr{Tx) {x = c oi VFT) 
and ar(Tx) (x = c or VFT) decay exponentially with 
d to a high approximation. The slower variation with 
T of Kr and and the overall values of Kr and ar all 
reflect decreased packing frustration with increasing d, 
thereby confirming the qualitative geometrical origin of 
the reduced fragility found upon increasing d. Both the 
reduced fragility with increasing d and the corresponding 
return to Arrhenius relaxation are consequences of the re¬ 
duction of packing frustration in high d where space is 
more “open” and packing constraints are more weakly 
“felt”. 


FIG. S4. (a) Main: Reduced thermal expansion coefficient 
ar — Tap as a function of T/Ta for various d. Crosses, 
asterisks, and circles iirdicate the points correspondiirg to Ta, 
Tc, and Tvft, respectively. Data are presented in the T range 
(Tvft, Ta) of interest in the present work. Inset: Same plot 
but shown in a linear scale representation, (b) Main: ariTA) 
as a function of d. Inset: ar{Tc) and ar{TvFT) as a function 
of d. Lines in the set indicate exponential fits: ar{Tr) = 
Ax exp(—Bajd) (x = c or VFT) with Ax and Bx the fitting 
parameters. We obtained the fitting estimates: Ac = 12.7 
and Be = 0.89 and Avft = 42.1 and Bvft = 1.63. 
















